Example 1 


The upward velocity of a rocket is given as a function of time in 
Table 1. 


Table 1 Velocity as a function of time. 


t(s) v(t) (m/s) 
0 0 

10 227.04 
15 362.78 
20 517.35 
22.5 602.97 
30 901.67 


Determine the distance, s, covered by the rocket from t =11 to 
t=16 using the velocity data provided and use any applicable 
numerical technique. 


Solution 


Method 1: Average Velocity Method 
The velocity of the rocket is not provided at t=11 and t =16, so 


we will have to use an interval that includes [11,16] to find the 
average velocity of the rocket within that range. In this case, the 
interval [10, 20] will suffice. 

v(10) = 227.04 

v(15) = 362.78 

v(20) = 517.35 
v(10) + v(15) + v(20) 

3 

_ 227.04 + 362.78 + 517.35 


3 
= 369.06 m/s 
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Figure 1 Velocity vs. time data for the rocket example 


Using 
s=VAt, 
we get 
s = (369.06)(16 —11) = 1845.3 m 


Method 2: Trapezoidal Rule 
If we were finding the distance traveled between times in the data 


table, we would simply find the area of the trapezoids with the 
corner points given as the velocity and time data points. For 


example 
20 


15 20 
f v(t)dt = f v(t)dt + f v(t)dt 
10 10 15 
and applying the trapezoidal rule over each of the above integrals 
gives 
20 
[v(ode = =— [v(10) + v(15)] + 
10 
The values of v(10), v(15) and v(20) are given in Table 1. 


However, we are interested in finding 
16 


15 16 
[ v(oat = f vedt + f (oat 
11 11 15 
and applying the trapezoidal rule over each of the above integrals 
gives 
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[v(15) + v(20)] 
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How do we find v(11) and v(16)? We use linear interpolation. 
To find v(11), 

v(t) = 227.04 + 27.148(t -10), 10 <t <15 

v(11) = 227.04 + 27.148(11 — 10) 

= 254.19 m/s 

and to find v(16) 

v(t) = 362.78 + 30.913(t -15), 15 < t < 20 

v(16) = 362.78 + 30.913(16 —15) 


= 393.69 m/s 
Then 
16 
15-11 16-15 
Í vdt = —— (11) + 362.78) + £ > (862.78 + v(16)) 
11 
_ 15-11 16-15 


(254.19 + 362.78) + 


(362.78 + 393.69) 
=1612.2 m 


Method 3: Polynomial interpolation to find the velocity profile 
Because we are finding the area under the curve from [10, 20} we 


must use three points, t =10, t=15, and t = 20, to fit a quadratic 


polynomial through the data. Using polynomial interpolation, our 
resulting velocity function is (refer to notes on direct method of 
interpolation) 


v(t) = 12.05 + 17.733t + 0.3766t°, 10 < t < 20. 


Now, we simply take the integral of the quadratic within our limits, 
giving us 


sx [2.05 +17.733t + 0.3766? hat 
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= 1604.3 m 


Method 4: Spline interpolation to find the velocity profile 

Fitting quadratic splines (refer to notes on spline method of 
interpolation) through the data results in the following set of 
quadratics 


v(t) = 22.704t, 0<t<10 
= 0.8888t" + 4.928t + 88.88, 10<t<15 
= —0.1356t° + 35.66t — 141.61, 15<t<20 
= 1.6048t? — 33.956t + 554.55, 20 <t < 22.5 
= 0.20889t?° + 28.86t — 152.13, 22.5<t<30 


The value of the integral would then simply be 


15 16 
s= [v(dt + | (oat 
11 15 


15 16 


~ f (0.8888? + 4.928t + 88.88)dt + f (-0.1356t? + 35.66 — 141.61 dt 
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141.61 
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= (16? -15?)-141.61(16-15) 


= 1595.9 m 


